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We study a problem of distributed detection of a stationary point event in a large extent wireless 
sensor network (WSN), where the event influences the observations of the sensors only in the vicinity of 
where it occurs. An event occurs at an unknown time and at a random location in the coverage region 
(or region of interest (ROI)) of the WSN. We consider a general sensing model in which the effect of the 
event at a sensor node depends on the distance between the event and the sensor node; in particular, in 
the Boolean sensing model, all sensors in a disk of a given radius around the event are equally affected. 
Following the prior work reported in |IT], ||2], [|3], the problem is formulated as that of detecting the event 

\o ; 

' and locating it to a subregion of the ROI as early as possible under the constraints that the average 

. run length to false alarm ('ARL2FAj is bounded below by 7, and the probability of false isolation (Pf\) 

' is bounded above by a, where 7 and a are target performance requirements. In this setting, we propose 

distributed procedures for event detection and isolation (namely MAX, ALL, and HALL), based on the 
local fusion of CUSUMs at the sensors. For these procedures, we obtain bounds on the maximum mean 
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detection/isolation delay (SADD), and on ARL2FA and PFI, and thus provide an upper bound on SADD 
as min{7, 1/a} — > oo. For the Boolean sensing model, we show that an asymptotic upper bound on the 
maximum mean detection/isolation delay of our distributed procedure scales with 7 and a in the same 
way as the asymptotically optimal centralised procedure ||2|. 

Index Terms 

Disorder problem, distributed quickest change detection, detection with distance dependent sensing, 
fusion of CUSUMs, multi-decision change-point detection, multi-hypothesis change detection 

I. Introduction 

Event detection is an important application for which a wireless sensor network (WSN) is deployed. A 
number of sensor nodes (or "motes") that can sense, compute, and communicate are deployed in a region 
of interest (ROI) in which the occurrence of an event (e.g., crack in a structure) has to be detected. In our 
work, we view an event as being associated with a change in the distribution ( or cumulative distribution 
function) of a physical quantity that is sensed by the sensor nodes. Thus, the work we present in this 
paper is in the framework of quickest detection of change in a random process. In the case of small 
extent networks, where the coverage of every sensor spans the whole ROI, and where we assume that 
an event affects all the sensor nodes in a statistically equivalent manner, we obtain the classical change 
detection problem whose solution is well known (see, for example, IH, Q, ||6l). In Q and lH, we have 
studied variations of the classical problem in the WSN context, where there is a wireless communication 
network between the sensors and the fusion centre lEl, and where there is a cost for taking sensor 
measurements lUl. 

However, in the case of large extent networks, where the ROI is large compared to the coverage region 
of a sensor, an event (e.g., a crack in a huge structure, gas leakage from a joint in a storage tank) affects 
sensors that are in its proximity; further the effect depends on the distances of the sensor nodes from 
the event. Since the location of the event is unknown, the post-change distribution of the observations 
of the sensor nodes are not known. In this paper, we are interested in obtaining procedures for detecting 
and locating an event in a large extent network. This problem is also referred to as change detection and 
isolation (see jU, ||2l, Q, Q, |[TOl ). Since the ROI is large, a large number of sensors are deployed to 
cover the ROI, making a centralised solution infeasible. In our work, we seek distributed algorithms for 
detecting and locating an event, with small detection delay, subject to constraints on false alarm and 
false isolation. The distributed algorithms require only local information from the neighborhood of each 
node. 
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A. Discussion of Related Literature 

The problem of sequential change detection/isolation with a finite set of post-change hypotheses was 
introduced by Nikiforov [1], where he studied the change detection/isolation problem with the observations 
being conditionally independent, and proposed a non-Bayesian procedure which is shown to be maximum 
mean detection/isolation delay optimal, as the average run lengths to false alarm and false isolation go 
to oo. Lai ifTOll considered the multi-hypothesis change detection/isolation problem with stationary pre- 
change and post-change observations, and obtained asymptotic lower bounds for the maximum mean 
detection/isolation delay. 

Nikiforov also studied a change detection/isolation problem under the average run length to false alarm 
(ARL2FA) and the probability of false isolation (PFI) constraints ||2], in which he showed that a CUSUM- 
like recursive procedure is asymptotically maximum mean detection/isolation delay optimal among the 
procedures that satisfy ARL2FA ^ 7 and PFI ^ a asymptotically, as min{7, 1/a} ^00 . Tartakovsky 
in 131 also studied the change detection/isolation problem where he proposed recursive matrix CUSUM 
and recursive matrix Shiryayev-Roberts tests, and showed that they are asymptotically maximum mean 
delay optimal over the constraints ARL2FA ^ 7 and PFI ^ a asymptotically, as min{7, 1/a} — > cxo. 

Malladi and Speyer [11] studied a Bayesian change detection/isolation problem and obtained a mean 
delay optimal centralised procedure which is a threshold based rule on the a posteriori probability of 
change corresponding to each post-change hypothesis. 

Centralised procedures incur high communication costs and distributed procedures would be desirable. 
In this paper, we study distributed procedures based on CUSUM detectors at the sensor nodes where the 
CUSUM detector at sensor node s is driven only by the observations made at node s. Also, in the case 
of large extent networks, the post-change distribution of the observations of a sensor node, in general, 
depends on the distance between the event and the sensor node which is unknown. 

B. Summary of Contributions 

1) As the WSN considered is of large extent, the post-change distribution is unknown, and could 
belong to a set of alternate hypotheses. In Section Jill we formulate the event detection/isolation 
problem in a large extent network in the framework of |0, ||3] as a maximum mean detec- 
tion/isolation delay minimisation problem subject to an average run length to false alarm (ARL2FA) 
and probability of false isolation (PFI) constraints. 

2) We propose distributed detection/isolation procedures MAX, ALL, and HALL (Hysteresis modified 
ALL) for large extent networks in Section |IV] The procedures MAX and ALL are extensions of 
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the decentralised procedures MAX 161 and ALL ||9l, |[T2l . which were developed for small extent 
networks. The distributed procedures are energy-efficient compared to the centralised procedures. 
Also, the known centralised procedures are applicable only for the Boolean sensing model. 
3) In Section ITVl we first obtain bounds on ARL2FA, PFL and maximum mean detection/isolation 
delay (SADD) for the distributed procedures MAX, ALL, and HALL. These bounds are then applied 
to get an upper bound on the SADD for the procedures when ARL2FA ^ 7, and PFI ^ a, 
where 7 and a are some performance requirements. For the case of the Boolean sensing model, 
we compare the SADD of the distributed procedures with that of Nikiforov's procedure |21 (a 
centralised asymptotically optimal procedure) and show that the an asymptotic upper bound on the 
maximum mean detection/isolation delay of our distributed procedure scales with 7 and a in the 
same way as that of 121. 

II. System Model 

Let ^ C be the region of interest (ROI) in which n sensor nodes are deployed. All nodes are 
equipped with the same type of sensor (e.g., acoustic). Let l^^'> G ^ be the location of sensor node s, and 
define £ := • • • We consider a discrete-time system, with the basic unit of time being 

one slot, indexed by A; = 0, 1, 2, • • • , the slot k being the time interval [k, k + 1). The sensor nodes are 
assumed to be time-synchronised (see, for example, |[T3l ). and at the beginning of every slot k^ \, each 

is) 

sensor node s samples its environment and obtains the observation Xj, G M. 
A. Change/Event Model 

An event (or change) occurs at an unknown time T G {1, 2, • • • } and at an unknown location G -A.. 
We consider only stationary (and permanent or persistent) point events, i.e., an event occurs at a point 
in the region of interest, and having occurred, stays there forever. Examples that would motivate such 
a model are 1) gas leakage in the wall of a large storage tank, 2) excessive strain at a point in a large 
2-dimensional structure. In |[T4l and lITSl . the authors study change detection problems in which the 
event stays only for a finite random amount of time. 

An event is viewed as a source of some physical signal that can be sensed by the sensor nodes. Let he 
be the signal strength of the evenj^. A sensor at a distance d from the event senses a signal hf.p{d) + W , 

*In case, the signal strength of the event is not laiown, but is known to lie in an interval \h, h\, we work with he — h as, 
this corresponds to the least Kullback-Leibler divergence between the "event not occurred" hypothesis and the ''event occurred' 
hypothesis. See 1 16] for change detection with unknown parameters for a collocated network. 
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where is a random zero mean noise, and p{d) is the distance dependent loss in signal strength which 
is a decreasing function of the distance d, with p{0) = 1. We assume an isotropic distance dependent 
loss model, whereby the signal received by all sensors at a distance d (from the event) is the same. 

Example 1 The Boolean model (see |17|): In this model, the signal strength that a sensor receives is 
the same (which is given by he) when the event occurs within a distance of from the sensor and is 
otherwise. Thus, for a Boolean sensing model. 



for some path loss exponent 77 > 0. For free space, ry = 2. 

B. Detection Region and Detection Partition 

In Example 2, we see that the signal from an event varies continuously over the region. Hence, unlike 
the Boolean model, there is no clear demarcation between the sensors that observe the event and those that 
do not. Thus, in order to facilitate the design of a distributed detection scheme with some performance 
guarantees, in the remainder of this section, we will define certain regions around each sensor. 

Definition 1 Given < /ii ^ he, the Detection Range r^i of a sensor is defined as the distance from 
the sensor within which the occurrence of an event induces a signal level of at least pi, i.e., 

rd := sup {d : hep{d) > pi} . 

■ 

In the above definition, pi is a design parameter that defines the acceptable detection delay. For a given 
signal strength he, a large value of pi results in a small detection range (as p{d) is non-increasing in 
d). We will see in Section HV-FI (Eqn. (fTTl) ) that the SADD of the distributed change detection/isolation 
procedures we propose, depends on the detection range r^, and that a small (i.e., a large ^1) results 
in a small SADD, while requiring more sensors to be deployed in order to achieve coverage of the ROI. 

We say that a location x E ROI is detection-covered by sensor node s, if p*^*^ — x|| ^ r^. For any 
sensor node s, V^'^^ := {x €^ A : — x|| ^ r^} is called its detection-coverage region (see Fig. [T]). We 
assume that the sensor deployment is such that every x ^ A is detection-covered by at least one sensor 



p{d) 




Example 2 The power law path-loss model (see |17|) is given by 



p{d) 



January 26, 2013 



DRAFT 



Fig. 1. Partitioning of ^ in a large WSN by detection regions: (a simple example) The coloured solid circles around each 
sensor node denote their detection regions. The four sensor nodes divide the ROI, indicated by the square region, into regions 
Ai, - ■ ■ , Ae such that region Ai is detection-covered by a unique set of sensors A/i. For example, Ai is detection covered by 
the set of sensors A/i = {1, 2, 4}, etc. 

(Fig. [Til. For each x £ A, define J\f{x) to be the largest set of sensors by which x is detection-covered, 
i.e., M{x) := {s : X e P^'*)}. Let C{M) = {J\f{x) : x £ A}. C{M) is a finite set and can have at most 
2" — 1 elements. Let N = \C{N')\. For each Mi G C{J\f), we denote the corresponding detection-covered 
region by Ai = A{Mi) := {x G ROI : 7V(x) = Mi}. Evidently, the Ai,l ^ i ^ partition the ROI. 
We say that the ROI is detection— partitioned into a minimum number of subregions, Ai,A2, ■ ■ ■ ,An, 
such that the subregion Ai is detection-covered by a unique set of sensors Mi, and Ai is the maximal 
detection-covered region of Mi, i.e., Vi / i'. Mi ^ Mi' and Ai fl Ai' = 0. See Fig. [T]for an example. 

C. Sensor Measurement Model 

is) 

Before change, i.e., for k < T, the obsei^vation at the sensor s is just the zero mean sensor 
noise VF^*\ the probability density function (pdf) of which is denoted by /o(-) {pre-change pdf). After 
change, i.e., for ^ T with the location of the event being i^, the observation of sensor s is given by 
xjf '' = hep{de,s) + wlf ^ where de,s '■= W^^^'^ —^e\\, the pdf of which is denoted by de,s) {post-change 
pdf). The noise processes {VF^^^} are independent and identically distributed (iid) across time and across 
sensor nodes. In the rest of the paper, we consider /o(-) to be Gaussian with mean and variance a^. 

We denote the probability measure when the change happens at time T and at location £e by P^*^^"^^ {•}, 
where d(^e) = [de,i: de,2, • • • ) de^n], and the corresponding expectation operator by E^^^"^'' [•]. In the case 
of Boolean sensing model, the post-change pdfs depend only on the detection subregion where the event 
occurs, and hence, we denote the probability measure when the event occurs at ie G Ai and at time T 
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by {■}, and the corresponding expectation operator by [•]. 
D. Local Change Detectors 

(s) 

We compute a CUSUM statistic C^, , /c ^ 1 at each sensor s based only on its own observations. The 
CUSUM procedure was proposed by Page [5| as a solution to the classical change detection problem 
(CDP, in which there is one pre-change hypothesis and only one post-change hypothesis). The optimality 
of CUSUM was shown for conditionally iid observations by Moustakides in ifTSl for a maximum mean 
delay metric introduced by PoUak 1191 which is SADD(t) := sup Et [t - T\t ^ T]. 

The driving term of CUSUM should be the log likelihood-ratio (LLR) of xjf^ defined as zjf\de^s) '■= 
('e)"' . As the location of the event £e is unknown, the distance de s is also unknown. Hence, 
one cannot work with the pdfs d^^s)- We propose to drive the CUSUM at each node s with Z^\rii), 

(s) 

where we recall that is the detection range of a sensor Based on the CUSUM statistic % A; ^ 1, 
sensor s computes a sequence of local decisions D^^^ S {0, 1}, A; ^ 1, where represents no-change 
and 1 represents change. For each set of sensor nodes Mi that detection partitions the ROI, we define 
^(M)^ the stopping time (based on the sequence of local decisions D^'^s for all s E Mi) at which the set 
of sensors Mi detects the event. The way we obtain the local decisions D^^^ from the CUSUM statistic 
Clf\k-^ 1, and the way these local decisions determine the stopping times t^-^'^\ varies from rule to 
rule. Specific rules for local decision and the fusion of local decisions will be described in Section JV] 
(also see [20]). 

An implementation strategy for our distributed event detection/isolation procedure can be the following. 
We assume that the sensors know to which detection sensor sets A/is they belong. This could be done 
by initial configuration or by self-organisation. When the local decision of sensor s is 1, it broadcasts 
this fact to all sensors in its detection neighbourhood. In practise, the broadcast range of these radios is 
substantially larger than the detection range. Hence, the local decision of s is learnt by all sensors s' that 
belong to Mi to which s belongs. When any node learns that all the sensors in Mi have reached the local 
decision 1, it transmits an alarm message to the base station lIlTI . A distributed leader election algorithm 
can be implemented so that only one, or a controlled number of alarms is sent. This alarm message is 
carried by geographical forwarding |22|- A system that utilises such local fusion (but with a different 
sensing and detection model) was developed by us and is reported in ll23l . 
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E. Influence Region 

After a set of nodes A/i declares an event, the event is isolated to a region associated with A/i called the 
influence region. In the Boolean sensing model, if an event occurs in Ai, then only the sensors s G Mi 
observe the event, while the other sensors s' ^ A/i only observe noise. On the other hand, in the power 
law path-loss model, sensors s' ^ A/i can also observe the event, and the driving term of the CUSUMs 
of sensors s' may be affected by the event. The mean of the driving term of CUSUM of any sensor s is 
given by 

Thus, the mean of the increment that drives CUSUM of node s decreases with d^^s and becomes negative 
when 2p{de,s) < p{fd)- In this region, we are interested in finding Te, the expected time for the CUSUM 
statistic Cl'^ to cross the threshold c. Define t^'^ := inf |A; : cjf^ ^ c}, and hence, Te = e['^^^'^'^ [t^'^ . 

Lemma 1 If the distance between sensor node s and the event, de^s i^ such that 2p{de,s) < pifd), then 

Te ^ exp(a;oc) 



where wq = 1 — 



Proof: From (Eqn. 5.2.79 pg. 177 of) |24|, we can show that E^'^*-^"'*'' [t^'^^] ^ exp(a;oc) where ujq 
is the solution to the equation 



0, 



p(d(4)) 

which is given hy ojq = \ — (see Eqn. ([1])). ■ 

We would be interested in Te ^ exp(a;o • c) for some < luq < 1. We now define the influence range 
of a sensor as follows. 

Definition 2 Influence Range of a sensor, r^, is defined as the distance from the sensor within which 
the occurrence of an event can be detected within a mean delay of exp (wgc) where loq is a parameter of 
interest and c is the threshold of the local CUSUM detector. Using Lemma [T] we see that = m.m{d' : 

2/o(d') ^ (I-^o)pM}- ■ 

A location x G ^ is influence covered by a sensor s if ll^*-*-* — x|| ^ r^, and a set of sensors J\fj is said 
to influence cover x if each sensor s G Mj influence covers x. 

From Lemma [H we see that by having a large value of loq, i.e., luq close to 1, the sensors that are 
beyond a distance of from the event take a long time to cross the threshold. However, we see from 
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(a) Detection and influence regions (b) Detection and influence regions 

of the Boolean model of the power law path loss model 

Fig. 2. Influence and detection regions: A simple example of partitioning of ^ in a large WSN. The coloured solid circles 
around each sensor node denote their detection regions. The four sensor nodes, in the figure, divide the ROI, indicated by the 
square region, into regions ^i, • ■ • , Ae such that region Ai is detection-covered by a unique set of sensors Mi. The dashed 
circles represent the influence regions. In the Boolean model, the influence region of a sensor coincides with its detection region. 



the definition of influence range that a large value of Wg gives a large influence range r^. We will see 
from the discussion in Section III-FI that a large influence range results in the isolation of the event to a 
large subregion of A. On the other hand, from Section IIV-EI we will see that a large loq decreases the 
probability of false isolation, a performance metric of change detection/isolation procedure, which we 
define in Section |llll 

We define the influence-region of sensor s as — x|| ^ Vi}. For the Boolean 

sensing model, = r^, and hence, V^^^ = T^^^ for all 1 ^ s ^ n, and for the power law path-loss 
sensing model, rj > r^, and hence, T>^^^ C T*^*-* for all 1 ^ s ^ n (see Fig. |2]). 

Recalling the sets of sensors Mi, 1 ^ i ^ N , defined in Section Hl-Bi we define the influence region 
of the set of sensors Mi as the region Bi such that each x £ Bi is within the influence range of all the 
sensors in A/",, i.e., B, := B{Mi) := RseM'^^'^- ^ote that A{Mi) = ( H ^^'M H ( (1 
where V is the complement of the set V, and C 7'^''). Hence, A{Mi) C B{Mi). For the power law 
path-loss sensing model, V^""^ C T''"^ for all 1 ^ s ^ n, and hence, A{Mi) C B{Mi) for all 1 5^ i ^ iV. 
For the Boolean sensing model, A{Mi) = B{Mi) fl ( fl V^^') ) , and hence A{Mi) = B{Mi) only when 
Mi = {1,2, ••• ,n}. Thus, for a general sensing model, A{Mi) C B{Mi). We note here that in the 
Boolean and the power law path loss models, an event which does not lie in the detection subregion of 
Mi, but lies in its influence subregion (i.e., £(, G B{Mi) \ A{Mi)) can be detected due to Mi because of 
the stochastic nature of the observations; in the power law path loss sensing model, this is also because 
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of the difference in losses p{de,s) between different sensors. 

Remark: The definition of the detection and influence ranges have involved two design parameters /xi 
and which can be used to "tune" the performance of the distributed detection schemes that we develop. 

■ 

F. Isolating the Event 

In Section II D, we provided an outline of a class of distributed detection procedures that will yield 
a stopping rule. On stopping, a decision for the location of the event is made, which is called isolation. 
In Section ITVl we will provide specific distributed detection/isolation procedures in which stopping will 
be due to one of the sensor sets Mi. 

An event occurring at location G -^i can influence sensors s' which influence cover i^, and hence, 
the detection can be due to sensors Mi ^ Mj which influence cover if.. Thus, we isolate the event to the 
influence region of the sensors that detect the event. Because of noise, detection can be due to a sensor 
set Mh which does not influence cover the event. Such an error event is called false isolation. 

An event occurring at £e G A is influence covered by sensors s' G M{£e) ■= {s : \\£^^^ — £e\\ ^ fi}. 
Hence, the detection due to any Mj Q M{ie) corresponds to the isolation of the event, and that due to 
Mj 2 M{£e) corresponds to false isolation. Note that in the case of Boolean sensing model M{£e) = Mi. 

In Section |llll we formulate the problem of quickest detection of an event and isolating the event to 
one of the influence subregions 81,82, - ■ ■ under a false alarm and false isolation constraint. 

III. Problem Formulation 

We are interested in studying the problem of distributed event detection/isolation in the setting devel- 
oped in Section |lll Given a sample node deployment (i.e., given £), and having chosen a value of the 
detection range, r^, we partition the ROI, A into the detection-subregions, Ai,A2, ■ ■ ■ ,An. Let Mi be 
the set of sensors that detection-cover the region Ai. Having chosen the influence range r^, the influence 
region Bi of the set of sensor nodes Mi can be obtained. We define the following set of hypotheses 

Ho : event not occurred, 

tlT,i '■ event occurred at time T in subregion Ai, T = 1,2, ■ ■ ■ , i = 1,2, - ■ ■ ,N. 

The event occurs in one of the detection subregions A-i, but we will only be able to isolate it to 
one of the influence subregions Bi that is consistent with the Ai (see Section III-FI) . We study distributed 
procedures described by a stopping time r, and an isolation decision L{t) € {1, 2, • • • , N} (i.e., the tuple 
{t,L)) that detect an event at time r and locate it to L{t) (i.e., to the influence region B^^^) subject 
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to a false alarm and false isolation constraint. The false alarm constraint considered is the average run 
length to false alarm ARL2FA, and the false isolation constraint considered is the probability of false 
isolation PFI, each of which we define as follows. 

Definition 3 The Average Run Lengtli to False Alarm ARL2FA of a change detection/isolation proce- 
dure r is defined as the expected number of samples taken under null hypothesis Hq to raise an alarm, 
i.e., 

ARL2FA(r) := E^, [r] , 

where Eoo[ ] is the expectation operator (with the corresponding probability measure being Poo{}) when 
the change occurs at infinity. ■ 

Definition 4 The Probability of False Isolation PFI of a change detection/isolation procedure r is 
defined as the supremum of the probabilities of making an incorrect isolation decision, i.e., 



>(d(4)) 
1 

where we recall that A/'(4) = {s : - 4|| ^ r^} is the set of sensors that influence covers ^ -^i- 



PFI(r) := max sup max P, iLir) = i} 



In the case of Boolean sensing model, the post-change pdfs depend only on the index i of the detection 
subregion where the event occurs, and hence, the PFI is given by 

PFI(r) := max max {L(t) = j} . 

In 121, Nikiforov defined the probability of false isolation, also, over the set of all possible change times, 
as SPFI(r) := supi^j^^y ^'^Pi^j^iiS.N sup^^^ P^^ {L{t) = j \ t ^ T}. Define the following classes 
of change detection/isolation procedures, 

A(7,a) := {(r, L) : ARL2FA(t) ^ 7, SPFI(t) ^ a} , 
A(7,a) := {(r, L) : ARL2FA(t) ^ 7, PFI(t) ^ a} . 

We define the supremum average detection delay SADD performance for the procedure r, in the same 
sense as PoUak |[T9l (also see 13), as the maximum mean number of samples taken under any hypothesis 
HT,i, ^ = 1) 2, • • • , A^, to raise an alarm, i.e., 

SADD(r) := sup sup E^.^^^^^^^ [r - T|r ^ T] . 
i.&A T^l 



January 26, 2013 



DRAFT 



12 



We are interested in obtaining an optimal procedure r that minimises the SADD subject to the average 
run length to false alarm and the probability of false isolation constraints, 

inf sup sup e5?^^=^^ [r - T|r ^ T] 

L€A T^l 

ARL2FA(r) ^ 7 

subject to 

PFI(r) ^ a. 

The change detection/isolation problem that we pose here is motivated by the framework of |[T], |I2], 
|[3]| . which we discuss in the next subsection. 

A. Centralised Recursive Solution for the Boolean Sensing Model 

In 121, Nikiforov and in |3|, Tartakovsky studied a change detection/isolation problem that involves 
> 1 post-change hypotheses (and one pre-change hypothesis). Thus, their formulation can be applied 
to our problem. But, in their model, the pdf of Xjt for k ^ T, under hypothesis gi is completely 

known. It should be noted that in our problem, in the case of power law path-loss sensing model, the 
pdf of the observations under any post-change hypothesis is unknown as the location of the event is 
unknown. The problem posed by Nikiforov [2| is 



r(0 

lili OUp OUp L 

■,L)eA(7,a) l^i<:N 

and that by Tartakovsky ||3l is 



inf sup sup E^^ [r - r|r ^ T] , (2) 

(T,L)eA(7,a) l^i<:N 



inf sup sup E^-* [r — T|r ^ T] . (3) 

(T,L)GA(7,a) l^i<:N 

Nikiforov ||2l and Tartakovsky iQ obtained asymptotically optimal centralised change detection/isolation 
procedures as min{7, ^} 00, the SADD of which is given by the following theorem. 

Theorem 1 (Nikiforov 03) For the N— hypotheses change detection/isolation problem (for the Boolean 
sensing model) defined in Eqn. the asymptotically maximum mean delay optimal detection/isolation 
procedure t* has the property, 



SADD(r*) ~ max<( 



In 7 — ln(a) I f 1 

as mm < 7, — > ^ 00, 



min KL(qi,qn)^ min KLioj^q^) 



a 



where KL{-,-) is the Kullback-Leibler divergence function, and gi is the pdf of the observation 'K/. for 
k ^ T under hypothesis 'H.T,i- H 

Remark: Since, A(7, a) C A(7, q), the asymptotic upper bound on SADD for r* is also an upper bound 
for the SADD over the set of procedures in A(7,a). 
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In the case of Boolean sensing model, for any post-change hypothesis only the set of sensor 

nodes that detection cover (which is the same as influence cover) the subregion Ai switch to a post- 
change pdf fi (and the distribution of other sensor nodes continues to be /g). Since the pdf of the 
sensor observations are conditionally i.i.d., the pdf of the observation vector, in the Boolean sensing 
model, corresponds to the post-change pdf gi of the centralised problem studied by Nikiforov [2 1 and by 
Tartakovsky [3 |. Thus, their problem directly applies to our setting with the Boolean sensing model. In 
our work, however, we propose algorithms for the change detection/isolation problem for the power law 
sensing model as well. Also, the procedures proposed by Nikiforov and by Tartakovsky are (while being 
recursive) centralised, whereas we propose distributed procedures which are computationally simple. 

In Section|IVl we propose distributed detection/isolation procedures MAX, HALL and ALL and analyse 
their false alarm (ARL2FA), false isolation (PFI) and the detection delay (SADD) properties. 

IV. Distributed Change Detection/Isolation Procedures 

In this section, we study the procedures MAX and ALL for change detection/isolation in a distributed 
setting. Also, we propose a distributed detection procedure "HALL," and analyse the SADD, the ARL2FA, 
and the PFI performance. 

A. The MAX Procedure 

Tartakovsky and Veeravalli proposed a decentralised procedure MAX for a collocated scenario in 
|[6ll . We extend the MAX procedure to a large WSN under the ARL2FA and PFI constraints. Recalling 
Section ini each sensor node i employs CUSUM for local change detection between pdfs /o and fi{-;rd). 
Let T^*) be the random time at which the CUSUM statistic of sensor node i crosses the threshold c for 

(i) 

the first time. At each time k, the local decision of sensor node i, Dj^ is defined as 



The global decision rule r declares an alarm at the earliest time slot k at which all sensor nodes 
j E Mi for some i = 1,2, - ■ ■ ,N have crossed the threshold c. Thus, 



i.e., the MAX procedure declares an alarm at the earliest time instant when the CUSUM statistic of all 
the sensor nodes Mi corresponding to hypothesis H^ j of some i have crossed the threshold at least once. 




1, for k ^ rW. 



for A; < T 



r 
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The isolation rule is L{t) = arg mini^j^Tvl''"'^''^^'^''^^}. i-e-> to declare that the event has occurred in the 
influence region -Bj^(r) = B{J^l{t)) corresponding to the set of sensors ■N'l{t) that raised the alarm. 

B. ALL Procedure 

Mei, 191, and Tartakovsky and Kim, 1251 . proposed a decentralised procedure ALL, again for a collocated 
network. We extend the ALL procedure to a large extent network under the ARL2FA and the PFI 
constraints. Here, each sensor node i employs CUSUM for local change detection between pdfs /o 
and fi{-;rii). Let cj^^ be the CUSUM statistic of sensor node i at time k. The CUSUM in the sensor 
nodes is allowed to run freely even after crossing the threshold c. Here, the local decision of sensor node 
i is 



1, if C'i'^ ^ c. 



The global decision rule r*'-'- declares an alarm at the earliest time slot k at which the local decision of 
all the sensor nodes corresponding to a set Mi, for some z = 1, 2, • • • , A^, are 1, i.e., 

^ALL,(M) {fc : D^^ = 1, Vj G AA^} = inf [k : c\^^ ^ c, Vj G 

r^LL min{rALL,(M)^i^^^^|, 

The isolation rule is L(r) = arg mini^j<jAr{r^'-'-'(-'^)}, i.e., to declare that the event has occurred in 
the influence region i3i(,-) = B{Ml[t-)) corresponding to the set of sensors N^r) that raised the alarm. 

C. HALL Procedure 

Motivated by ALL, and the fact that sensor noise can make the CUSUM statistic fluctuate around the 
threshold, we propose a local decision rule which is when the CUSUM statistic has visited zero and 
has not crossed the threshold yet and is 1 otherwise. We explain the HALL procedure below. 
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(0 



The following discussion is illustrated in Fig. |3] Each sensor node i computes a CUSUM statistic Cj 
based on the LLR of its own observations between the pdfs fi{-;rd) and /q. Define Uq^ := 0. Define 
as the time at which Cj. crosses the threshold c (for the first time) as: 

V^^ :=inf |A: : C^^ ^ c} 

(see Fig. |3] where the "overshoots" R^^\ at V^^\ are also shown). Note that inf := oo. Next define 

[/« :=inf{A;>y/*) :Ci')=0}. 

Now starting with u\^\ we can recursively define \ etc. in the obvious manner (see Fig. |3]l. Each 
node i computes the local decision D^^^ based on the CUSUM statistic C^'^ as follows: 

f 1, if V^'^ s^k< Uf for some j 
of = } ' ^ J (4) 

I 0, otherwise. 

The global decision rule is a stopping time r^^^^ defined as the earliest time slot k at which all the 
sensor nodes in a region have a local decision 1, i.e., 

The isolation rule is L{t) = &igmSxii<^i<^N{T^^^^''^^''>}, i.e., to declare that the event has occurred in 
the influence region ^l(t) = '^(-^l(t)) corresponding to the set of sensors M^t) that raised the alarm. 
Remark: The procedures HALL, MAX and ALL differ only in their local decision rule; the global decision 
rule as a function of {-D^*''} is the same for HALL, MAX and ALL. For the distributed procedures MAX, 
ALL, and HALL, we analyse the ARL2FA in Section HV-Dl the PFI in Section ITV-El and the SADD 
performance in Section IIV-FI 

D. Average Run Length to False Alarm ('ARL2FA) 

From the previous sections, we see that the stopping time of any procedure (MAX, ALL, or HALL) is 
the minimum of the stopping times corresponding to each Mi, i.e., 

^procedure ._ ^^-^ | ^ proced u re, (M) ; 1 ^ ^ ^ ivj . 

Under the null hypothesis Hq, the CUSUM statistics C^^''^ of sensors s ^ Mi are driven by independent 

(s) 

noise processes, and hence, s are independent. But, there can be a sensor that is common to two 
different A/is, and hence, 7-P''ocedure,(M)s^ in general, are not independent. We provide asymptotic lower 
bounds for the ARL2FA for MAX, HALL, and ALL, in the following theorem. 
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Theorem 2 For local CUSUM threshold c, 

ARL2FA(r^^^) ^ exp (aMAXc) • (1 + o(l)) (5) 
ARL2FA(rHALL) ^ exp (aHALLc) • (1 + o(l)) (6) 
ARL2FA(r^'-L) ^ exp (aALLc) • (1 + o(l)) (7) 

(0(1) — 7- fl5 c — > ooj, where for any arbitrarily small (5 > 0, cmax = ^HALL = 1 — (5, oall = rn — 5, 

where m = niin{7Vi \ |Jj_^j jex-^j • ^ ^ -^i' '^'^^ °f indices of the detection sets that are minimal 
in the partially order of set inclusion among the detection sets. 

Proof: See Appendix lAl ■ 
Thus, for procedure, for a given ARL2FA requirement of 7, it is sufficient to choose the threshold c as 



In 7 



(1 + 0(1)), as7^cDo. (8) 



E. Probability of False Isolation (PFI) 

A false isolation occurs when the hypothesis Hj- j is true for some i and the hypothesis j 7^ H^,! 
is declared to be true at the time of alarm, and the event does not lie in the region B{Mj). The following 
theorem provide asymptotic upper bounds for the PFI for each of the procedures MAX, ALL, and HALL. 

Theorem 3 For local CUSUM threshold c, 



pPI(^MAX) ^ ^^Piz^M^ . (1 + ,(1)) (9) 
-DM AX 

PFI(rHALL) ^ !^P(z^y^ . (1 + ^(1)) (10) 
-Bhall 

pPK^ALL) ^ + (11) 

-DALL 

where o(l) [) as c ^ 00, and omax = ohall = -^f^ oall = — n' ~ 

Boolean sensing model, ^ is 2 for Boolean sensing model and is 1 for path-loss sensing model, m = 
Tii\n{\Nj\M{le)\ : 1 ^ i ^ iV,4 e A, 1 J iV,A/;- 2 AA(4)} and 

m = max{|AAj \ A/'(4)| : 1 ^ « ^ iV, 4 G A, 1 ^ J ^ N,Nj % M{le)\, and Bmax. ^hall, and 5all 
are positive constants. 

Proof: See Appendix iBl ■ 
Thus, for a given PFI requirement of a, the threshold c for should satisfy 
_ -lni?procedure - Ina ,,,, - In a 



(1 + 0(1)) = (1 + 0(1)), asa^O. (12) 



^procedure ^procedure 
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F. Supremum Average Detection Delay ( SADD) 

In this section, we analyse the SADD perfomiance of the distributed detection/isolation procedures. 
We observe that for any sample path of the observation process, for the same threshold c, the MAX rule 
raises an alarm first, followed by the HALL rule, and then by the ALL rule. This ordering is due to the 
following reason. For each sensor node s, let r'-'*^ be the first time instant at which the CUSUM statistic 
C^*^ crosses the threshold c (denoted by y/*^ in Figure |3]l. Before time r^*^ the local decision is for 
all the procedures, MAX, ALL, and HALL. For MAX, for all k ^ t'^^\ the local decision D^^'' = 1. 
Thus, the stopping time of MAX is at least as early as that of HALL and ALL. The local decision of 

is) is) is) 

ALL is 1 (Z)^ = 1) only at those times k for which ^ c. However, even when Cf, <c, the local 



decision of HALL is 1 if V^'^^ ^ k < U^^' (see Figure O for some j. Thus, the local decisions of MAX, 



HALL, and ALL are ordered as, for all /fc ^ 1, L>^''^(MAX) ^ L'['^(HALL) ^ Dl"'{Ml), and hence, 
^MAX,(M) ^ ^HALL,(AA,) ^ ^ALL,(Arj g^^j^ stopping times MAX, HALL, or ALL is the minimum 

of stopping times corresponding to the sets of sensors {Mi : i = 1,2, ■ ■ ■ , N}, i.e., 

^procedure ^ j^^- j^|^procedure,(M) : i = 1, 2, ■ ■ • , A^} 



where "procedure" can be MAX or HALL or ALL. Hence, we have 



_MAX ^ ^HALL ^ ALL 



(13) 



From iH, we see that 



sup E^' 

T>1 



■^ALL,(M) _ I ^ALL,{K) > jn' 



J (1+0(1)) 



(14) 



where / is the KuUback-Leibler divergence between the post-change and the pre-change pdfs. For 
ie G Ai, we have Vs G Mi, de s ^ rd- Also, since r'^'-'- 5^ t-ALL,(a/;)^ j^^^g 



sup sup 

4eA T^ji 



ALL 



T I r^"-"- > T 



^ sup sup 



;^ALL,(M) 



T I r^"-"- > T 



(15) 



From Appendix O Eqn. ([TSl l becomes, 

„ALL 



sup sup 

4eA T>i 



T I r^'-'- > T 



^ sup sup 

4gA T>1 



KL(/i(.;rd),/o) 
From the above equation, and from Eqn. ([T3] l. we have 

SADD(r^^^) ^ SADD(r^^'-L) ^ SADD(r^LL) ^ 



KL(/i(.;rrf),/o) 



fo(l)) (16) 
(1 + 0(1)), as c 00, (17) 
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Remark: Recall from Section III-BI that /ii = hep{rd)- We now see that /xi govems the detection 
delay perfonnance, and fii can be chosen such that a requirement on SADD is met. Thus, to achieve a 
requirement on SADD, we need to choose appropriately. A small value of rd (gives a large /zi and 
hence,) gives less detection delay compared to a large value of r^. But, a small requires more sensors 
to detection-cover the ROI. 

In the next subsection, we discuss the asymptotic minimax delay optimality of the distributed procedures 
in relation to Theorem [T] 



G. Asymptotic Upper Bound on SADD 

For any change detection/isolation procedure to achieve a ARL2FA requirement of 7 and PFI require- 
ment of a, a threshold c is chosen such that it satisfies Eqns. [8]and[T2l i.e., 

c = max I ^^^,-^^^1 (1 + 0(1)). (18) 



'2 procedure ^procedure 

Therefore, from Eqn.dTTTi. the SADD is given by 

SADD(rP--^"-) ^ .niax|^^,—i^| (1 + 0(1)). (19) 

'^'-(Ji(-, Trfj, /oj (Oprocedure "proved ure J 

where o(l) — )• as min{7, ^} — )• 00. Note that as decreases, KL(/i(-; r^), /o) = ^^2^2'^ increases. 
Thus, to achieve a smaller detection delay, the detection range can be decreased, and the number of 
sensors n can be increased to cover the ROI. 

We can compare the asymptotic SADD performance of the distributed procedures HALL, MAX and 
ALL against Theorem [T] for the Boolean sensing model. For Gaussian pdfs /o and /i, the KL divergence 
between the hypotheses H-p j and Htj is given by 



(^')) dx 



1- f n W^] -i^n n m^^-^) n /o(-^^'^) 

= J;KL(/i,/o)- '<L(/i,/o)+ KL(/i,/o) 

= \MiAMj\ KL(/i,/o) 

where the operator A represents the symmetric difference between the sets. Thus, from Theorem [T] for 
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(20) 



Gaussian /o and /i, we have 

where a* = min IMA, 

and b* = min iTViAA/", !. 

The SADD performance of the distributed procedure with the Boolean sensing model is 

SADD(^procedure) ^ ^ ^lax | — | (1 + o(l)). 

'^'-[JliJO) I, C.pi-ocedure "procedure J 

where o(l) — ?> as min{7, ^} — > cx). Thus, the asymptotically optimal upper bound on SADD (which 
corresponds to the optimum centraUsed procedure r*) and that of the distributed procedures ALL, HALL, 
and MAX scale in the same way as ln7/KL(/i, /q) and — lna/KL(/i, /q). 

V. Numerical Results 

We consider a deployment of 7 nodes with the detection range = 1, in a hexagonal ROI (see Fig. ^ 
such that we get = 12 detection subregions, and Mi = {1, 3, 4, 6}, M = {1, 3, 4}, 7V3 = {1, 2, 3, 4}, 
Mi = {1,2,4}, M5 = {1,2,4,5}, Me = {2,4,5}, M7 = {2,4,5,7}, Ms = {4,5,7}, Mg = {4,5,6,7}, 
Mio = {4, 6, 7}, Mil = {3, 4, 6, 7}, and Mu = {3, 4, 6}. The pre-change pdf considered is /o ~ M{0, 1), 
and the detection range and the influence range considered are = 1.0 and = 1.5 respectively. 

We compute the SADD, the ARL2FA and the PFI performance of MAX, HALL, ALL, and Nikiforov's 
procedure (||2]) for the Boolean sensing model with fi ~ M{1, 1), and plot the SADD vs log(ARL2FA) 
performance in Fig. 5(a)| of the change detection/isolation procedures for PFI ^ 5 x 10^^. The local 



CUSUM threshold c that yields the target ARL2FA and other simulation parameters and results are 
tabulated in Table Jl To obtain the SADD the event is assumed to occur at time 1, which corresponds 
to the maximum mean delay (see ifTOl . |[26l ). We observe from Fig. |5(a)| that the SADD performance 



of MAX is the worst and that of Nikiforov's is the best. Also, we note that the performance of the 
distributed procedures, ALL and HALL, are very close to that of the optimal centralised procedure. 
For eg., for a requirement of ARL2FA = 10^ (and PFI ^ 5 x 10^^), we observe from Fig. |5(a)| that 
SADD(t^^^) = 26.43, SADD(rHALL) ^ 13 73^ SADD(rA'-L) ^ ^2.20, and SADD(t*) = 11.28. Since 

(s) 

MAX does not make use of the the dynamics of beyond r*, it's SADD vs ARL2FA performance is 

(s) 

poor. On the other hand, ALL and HALL make use of for all k and hence, give a better performance. 

For the same sensor deployment in Fig. |4l we compute the SADD and the ARL2FA for the square law 
path loss (rj = 2) sensing model given in Section JIl Also, the signal strength is taken to be unity. 
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Fig. 4. Sensor nodes placement: 7 sensor nodes (which are numbered 1,2, - • • ,7) represented by small filled circles are placed 
in the hexagonal ROI A. The sensor nodes partition the ROI into the detection subregions Ai,A2,- ' ,Ai2 (for both the 
Boolean and the power law path loss sensing models). 



TABLE I 

Simulation parameters and results for the Boolean sensing model for PFI 5 x 10 



Detection/ 
Isolation 
procedure 


No. of 

MC 

runs 


Threshold 

c 


ARL2FA 


99% Confidence interval 


SADD 


99% Confidence interval 


ARL2FAtow<:r 


A R L2 FAuppei- 


SADDlower 


SAD Dupper 


MAX 


10* 


2.71 


10^ 


93.69 


106.61 


8.77 


8.45 


9.09 


10* 


4.93 


10^ 


942.10 


1065.81 


14.89 


14.41 


15.37 


10* 


7.24 


10* 


9398.61 


10640.99 


21.01 


20.42 


21.61 


10* 


9.52 


10^ 


95696.90 


108008.89 


26.43 


25.76 


27.11 


HALL 


10* 


1.67 


10^ 


92.67 


107.58 


5.96 


5.72 


6.20 


10* 


2.69 


10=^ 


927.17 


1085.48 


8.81 


8.48 


9.14 


10* 


3.66 


10* 


9239.97 


10826.71 


11.58 


11.17 


11.99 


10* 


4.52 


10^ 


92492.85 


108389.15 


13.78 


13.32 


14.23 


ALL 


10* 


2.16 


10^ 


915.94 


1089.33 


7.82 


7.53 


8.11 


10* 


2.96 


10* 


9197.23 


10811.90 


10.07 


9.70 


10.44 


10* 


3.71 


10^ 


92205.45 


107952.43 


12.20 


11.76 


12.63 


Nikiforov 


10* 


2.75 


10^ 


98.30 


116.32 


4.75 


4.52 


4.98 


10* 


4.50 


10^ 


986.48 


1048.23 


7.08 


6.79 


7.38 


10* 


6.32 


10* 


9727.19 


10261.94 


9.14 


9.00 


9.68 


10* 


8.32 


10' 


98961.41 


110415.50 


11.28 


11.00 


12.25 
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ARL2FA (samples) ARL2FA (samples) 

(a) SADD vs ARL2FA for the Boolean model (b) SADD vs ARL2FA for the square law path loss model 

Fig. 5. SADD versus ARL2FA (for PFI < 5 x 10"^) for MAX, HALL, ALL and Nikiforov's procedure for the Boolean and 

the square law path loss sensing models. In the Boolean sensing model, the system parameters are /o ~ A'^(0, 1), /i ~ A'^(0, 1), 
and in the case of path loss sensing model, the parameters are /o ~ A'^(0, 1), he = 1, Vd = 1.0, r-i = 1.5. 



TABLE II 

Simulation parameters and results for the square law path loss sensing model for PFI < 5 x 10 



Detection/ 
Isolation 
procedure 


No. of 

MC 

runs 


Threshold 

c 


ARL2FA 


99% Confidence interval 


SADD 


99% Confidence interval 


ARL2FAlo„er 


A R L2 FAypper 


SADDlower 


SAD Dupper 


MAX 


10* 


2.71 


10^ 


93.69 


106.61 


30.74 


29.31 


32.17 


10* 


4.93 


10^ 


942.10 


1065.81 


79.60 


75.86 


83.34 


10* 


7.23 


10* 


9398.61 


10640.99 


169.63 


161.61 


177.65 


10* 


9.52 


10^ 


95696.90 


108008.89 


301.77 


286.88 


316.66 


HALL 


10* 


1.67 


10^ 


92.67 


107.58 


20.58 


19.43 


21.74 


10* 


2.69 


10^ 


927.17 


1085.48 


40.56 


38.24 


42.88 


10* 


3.66 


10* 


9239.97 


10826.71 


66.45 


62.57 


70.33 


10* 


4.52 


10^ 


92492.85 


108389.15 


96.93 


91.03 


102.82 


ALL 


10* 


1.33 


10^ 


92.24 


107.79 


20.19 


19.06 


21.32 


10* 


2.16 


10^ 


915.94 


1089.33 


39.90 


37.59 


42.21 


10* 


2.96 


10* 


9197.23 


10811.90 


63.34 


59.43 


67.24 


10* 


3.71 


10' 


92205.45 


107952.43 


98.96 


93.01 


104.92 



Thus, the sensor sets (TVjS) and the detection subregions {AiS) are the same as in the Boolean model, we 
described above. Since is taken as 1, ~ -^(1; 1)- Thus, the LLR of observation X^^ is given 

by In ("^y^T^iy^) ~ -^k^ ~ \' which is the same as that in the Boolean sensing model. Hence, under the 
event not occurred hypothesis, the ARL2FA under the path loss sensing model is the same as that of the 
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Boolean sensing model. The CUSUM threshold c that yields the target ARL2FAs and other parameters 
and results are tabulated in Table |ll] To obtain the SADD the event is assumed to occur at time 1, and 
at a distance of rj from all the nodes of Mi that influence covers the event (which corresponds to the 



maximum detection delay). We plot the SADD vs log(ARL2FA) in Fig. 5(b) The ordering on SADD for 

any ARL2FA across the procedures is the same as that in the Boolean model, and can be explained in 

the same manner. The ambiguity in lf> affects /i(-;de,s) and shows up as large SADD values. 

VI. Conclusion 

We consider the quickest distributed event detection/isolation problem in a large extent WSN with a 
practical sensing model which incorporates the reduction in signal strength with distance. We formulate 
the change detection/isolation problem in the optimality framework of ||2] and ||3]. We propose distributed 
detection/isolation procedures, MAX, ALL and HALL and show that as min{ARL2FA, 1/PFI} — )• oo, the 
SADD performance of the distributed procedures grows in the same scale as that of the optimal centralised 
procedure of Tartakovsky [3 1 and Nikiforov 

Appendix A 
Proof of Theorem [2] 

From detection sensor sets A/i, i = 1, 2, • • • ,N, we choose the collection of indices X C {1, 2, ■ ■ • , N} 
such that any two sensor sets Mi, Mj, i, j ^ X, are not partially ordered by set inclusion. For each i ^ I, 
define the set of sensors that are unique to the sensor set Mi, Aii := Mi\ [j Mj C Mi. The sets 
Mli,A42, • • • ,Mt\x\ are disjoint. Under the null hypothesis, Hq, the observations of sensors in the sensor 
sets Aii,M-2, ■ ■ ■ ,Ml\x\ are iid, with the pdf /q ~ M{0,a'^). For every Mi, there exists Mj such that 
Mj C Mi, so that r™'^'(-^') ^ ^ru\e,(M,)_ Hence, rP''°'^^^^''^ = min{rP''°"^"''^'(-^-) : i = 1, 2, • • • , iV} ^ 

j^-j^|^procedure,(A^,) ; ^ g J} =; f rule_ ^l^^^e, 

ruiel ^ . p |^ rule ^ ^mcj jj^^ Markov inequality) 



-oo 



^rule 



T 



or, ^ ^ P |f ^"'^ > e™"| = n > ^"^ / • (21) 



We analyse Poo {r > e"^"} as c ^ oo, for ALL, MAX, and HALL. For ALL, 

p^|^ALL,(A^O = A;}^ Poo{c(^^^c,V.GA1,} = nP°o{4^^^c} 

^ g-cm, (using Wald's inequaUty) 
Therefore, Poo {r all,(A^O ^ A;} ^ A; • e"^"'' 



Poo {r ^'-'-.(A^O > e--} ^ 1 



g— c(mi— m) 
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Hence, for any m < mi, we have lim inf c_5>oo Poo {''" ^'-'-'(-^■) > e™^} = 1. A large m (which is smaller 
than all m-js) is desirable. Thus, a good choice for m is oall = minjmi : i G Z} — 5. for some arbitrarily 
small 5 > 0. Hence, from Eqn. (|2T]) . 



ALL 



^ exp (oallc) (1 + o(l)) 



(22) 



For MAX, at the stopping time of MAX, at least one of the CUSUM statistics is above the threshold c. 

Poo [r = A;} ^ Poo [C^^^ ^ c, for some s G M,] 



^mie ^ (using Wald's inequality). 



(23) 



Therefore, for any arbitrarily small 5 > 0, Poo \r '^^^■(^') > e^^""^)"^! ^ 1 - m^e 

liminfPoolT^^>^'(^')>e(i-'^)^| = l. 

c— ^oo I J 

Let omax = 1 — 5- For any arbitrarily small J > 0, we see from Eqn. (|2TI ). 



5c 



MAX 



^ exp((l -5)c) (l + o(l)) =: exp(aMAXc)(l + o(l)), 



(24) 



(25) 



For HALL, for the same threshold c, the stopping time of HALL is after that of MAX. Hence, r ^ 
r ^A^. Hence, Eoo [r ^^"-Lj ^ Eoo [r ^^^] ^ exp ((1 - 5)c) (1 + o(l)) (from Eqn. (gS])). Thus, for 
Call := 1 — 5, for any arbitrarily small 5 > 0, 



hall 



^ exp (aALLc) (1 + o(l)) 



(26) 



Appendix B 
Proof of Theorem [3] 

Consider l^. G Ai. The probability of false isolation when the detection is due to J\fj ^ N{le) is 

p(d{4)) l^rule ^ ^rule,(AG)| ^ p(d(4)) |^mle,(Ar,) ^ ^mle,(Ar^) ^ = 1,2,--- , ivj 

^ p(d(4)) |^mle,{Ar,) ^ ^mle,{M)| 

oo 

^ p{d{4)) |^rule,(M) ^ ^1 p(d(^=)) |^rule,(A/-,) ^ ^ | ^rule,(M) ^ ^| 

A: 



fc=l 

oo 



^ p{d{4)) |^rule,(M) = /c| ^ p(d(^=)) |^rule,(A/;) ^ ^ | ^rule,(M) ^ /.j 
fc=l Li=l 
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A. PFI(r'^'-'-) - Boolean Sensing Model 

^ exp {—\^fj \Mi\c) (using Wald's inequality). 



Therefore, jr^LL.CAG) ^ ^all,(M)| ^ exp{-\Mj \K\c) ■ E^^^ 

^exp{-{\Mj\K\c-Hc))) 



^ALL,(M) 
1 



■(l+o(l)). 



Hence, PFI(r'^'-L) ^ ^^^^^^ , . max P^') jr^L^CAG) ^ ^all,(a/;)| 



exp (— (mc — ln(c))) 



(27) 



where n = minjlA/'i : i = 1,2, ■ ■ ■ , N\}, m = min {lA/'j \ A/il}. For any n, there exists 

co(n) such that for all c> co{n),c < e*^/". Using this inequality, for sufficiently large c 



PFI(rALL)^!^(-((^-^)c)) 



na 



where 6all = m — 1/n and i?ALL = ZIck- 



B. PFI(r^^^) - Boolean Sensing Model 

pi^) |^MAX,(Ar,) ^ ^ I ^MAX,(M) = ^| ^ p« |^(.) ^ t,ys e Nj \ r^^^'^"^') = A:} 

^ Poo {r(") ^ i, Vs e Nj \Mi I t'^a^'(-^') = A;} 
t 

„MAx,(jv;) ^ ^1 



= n EPoo{rW = 

seMj\Mi n=i 



n T 



s€Afj\Mi n=l 



^ exp {-mjic) t"^'' 
Hence, P? {r^^^^'^-^) ^ r^^^^'^"^')} ^ exp(-m,,c) [(r'^^^^'-^-) 



^ exp (-mj,c) ^^yt;;^(1 + ^(l)) 
_ exp (- {rrijiC - (1 + mjj) ln(c))) 



(1 + «(!)) 
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Let m = mill mji, m = max and a* = min a^+"^-'* . 

Therefore, PFI(r^A^) ^ max max P^^^ /^max.ia/;) ^ ^max,(M)| 
^ exp (- (mc - (1 + m) ln(c))) ^ ^^^^^ 

For any n, there exists co(n) such that for all c > co(n), c < e'^/". Hence, for sufficiently large c 
PFI(t^a>^) ^ max max P« |r^A>^'(^) ^ r^f^M^^)] 

(1+0(1)) = 2^Et^^(i+o(i)). 

a -DMAX 

where 6max = m — ((1 + fh)/n) and 5max = «*• 

C. PFI(r"*'-'-) - Boolean Sensing Model 

pW |^HALL,(AG) = i I ^HALL,(M) = ^ pW |^(.) <^t,\/s^ Mj | r^^^'-^-^*) = k] 

which has the same form as that of MAX. Hence, from the analysis of MAX, it follows that 

p« |^HALL,(A^) ^ ^HALL,(M)| ^ exp (-m,,c) E« [(^HALL,(M))l+m,.- 



|7Vj|l+mj.Q,l+mj, 



^ exp i-rrijic) , + ^^(1)) 



= exp (- (mjjC - (1 + ruji) ln(c))) 



«|M|J 



(1 + «(!)) 



PFKr^^^LL) ^ p» r HALL,(A/:)^^HALL,(M)\ 

^ exp (- (mc - (1 + ?7t) ln(c))) 

where a* = min (a • \Afi\V~^"^'\ For any n there exists cnfn) such that for all c > 

co(n),c < e'^/". Hence, for sufficiently large c 

PFKr^^^LL) ^ p« f^HALL,(A/;)^^HALL,(M)j 

g "''<'<"T'^>'> (1+0(1)) = "'"-''"'^"•^> (1+0(1)). 

a -DHALL 

where 6hall = m — {{1 + m) /n) and -Bhall = «*• 



January 26, 2013 



DRAFT 



26 



PFI - Path-Loss Sensing Model 

Lemma 2 For s e J\fj \Me and for t ^ T, (with the pre-change pdf fo ~ M(0, a^) and the post-change 
pdfh^Af{hep{rs),a^)) 

w/zere we recall that the parameter defines the influence range, and a = KL{fi, /o). 
Proof: For seJ\fj\ Me and for t ^ T, 

p(d(^e)) |^(.) ^ 

= Pf'^f^^" I max y In f /ll^V^ll^ ^ c 



/c=l 



n = l 



^Ep4E...f^).f|:pr"{|:...(^).. 



n=l kfe=l 



n=i U=i \ /o(.Aj. j y J (^j.^-^ y /o(Aj, j y ^^j, 

T— 1 I'n 2 ^oc|'ri 2 / 



^ -/>(<(.,.) I +( 



^^Poo<{cxp(6i^X^"'j > exp ^^-^-^(c + nat^o)^ i for any (9 > 0. 



fe=i 

.2 



Hence, p(''(^^» {c}^' > c} ^ ^^exp (^--^(c + naa;j) (e 



2 

Since the above inequality holds for any > 0, we have 



The minimising 6 is Therefore, for 9 = 



n=l ^ 



— (c + nauio) 



Note that - -^^ , ~" ^ + = =^ + 



4Q;n 4a{n — 1) 4 4a(n — l)n 
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Therefore, by iteratively computing the exponent, we have 

< exp f-fcii^) ^ exp f-^^ln - 1)) exp ( l-) 

OO 



4a / V 4 / V 4a / 



L^^P z;;^ ^ «^p [-^^' 



n=l 



aoj' 



D. PFI(t^'-'-) - Path Loss Sensing Model 

pWi.)) |^ALL,(^) ^ ^ I ^ALL,(M) = ^1 ^ p(d(4)) ^ G AG- I rA^L,(M) = ^1 



^ p 



we.)) ^ ^^^^ g A/-, \ Ar(4) I T^'''(^') = A;} 

sGA/;\Ar(£e) 



|AGW(^e)| 



(from Lemma |2|) 



Therefore, P[''^'^^^ i^r^^UM) ^ ^all,(M)| ^ ^|A^AAr(£»)lE(d(^=)) 



^all,{m; 



a|7Vi[ 



Let m = mill lA/",- \ M(ie)\ and n = min||A/i| : i = 1,2, - ■ ■ , N}. Define K 

\mm^)\ 



max 



exp 



1— exp — 



. Therefore, 



PFI ( T^LL\ ^ p(d(^=)) r ALL,(^) ^ ^ALL,(M)1 



an 



-(1 + 0(1)). 

For any n there exists co(n) such that for all c > co(n), c < e'^^"'. Hence, for sufficiently large c 

PF,f,ALL^^^^fp(zSi^iMa + ofi)) = ^^py^tM:^(i + o(i)) 



an 



5ALL,d 



where ?)ALL,d = ("ia;o/2) - (1/n) and 5ALL,<i = an/K. 
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E. PFI(r^^^) - Path Loss Sensing Model 

p(d(4)) |^MAX,(A^) ^ ^ I ^MAX,(M) ^ ^1 <^ p(d(4)) ^ V s e AG- \M{Q I r^A^'(^') = A:} 

^ -Q p(d(4)) |^{.) ^ ^ [ ^MAX,(M) ^ ^1 

sGA/'AA^{4) 

sGArj\A/'(^») n=l 



^ . i|Ar,W(4)| (from LemmaHl) 



p(d(^e)) |^lvlH/^,^vvjj ^ 



.MAX,(A/^) ^ ^MAX,(M 



.MAX,(M)\l+|AOW(f.)| 



|A/'AA/'(4)| 



Let m 



min 



|AA,- \AA(4)|, m 



.i+|A/;-W(4)| 



a 



1+|AGW(4)| 



(1 + 0(1)) 



max 



define K = max 

.MAXn 



exp 



1— exp 



. Therefore, 



\M,\M{ie)\, and 



PFI(r'^^^)^ max sup max pW^O) f^MAX,(A/;) ^ ^MAX 



^ ^exp 
a* 



where a* 



mill qI+IA/'j Pqj- ^j^y ^ there exists co(n) such that for all 

c > co(n), c < e'^/". Hence, for sufficiently large c 



PFIfr 



MAX^ 



K 

^ ^exp 

a* 

rriLj,, 



.0 



1 + m 



where 6MAX,d = (^) - ( 



2 n 
and Bmaxa = 



c (l + o(l)) 



exp(-bMAX,d • c) ^ 

-BMAX,d 



l + o(l)), 



F. PFI(r^^^'-'-) - Path Loss Sensing Model 

p(d(4)) |^HALL,(A/-) ^ ^ I ^HALL,(M) ^ fcj ^ p(d(4)) |^(.) ^ ty , ^ J^^ | ^HALL,(M) ^ fcj 

which has the same form as that of MAX. Hence, from the analysis of MAX, it follows that 



5(d(4)) |, 



.HALL, (A/]) ^ ^HALL,(A/') 



Ar,W(4)lF(d(^'')) 



HALL,(M)\1+|A(,W(4)I 



l+|Ar,-W(4)| 



Therefore, PFIfr^^'-'-) ^ max sup max p(d(^0) f. 



.HALL,(A/^) ^ .j.HALL,(M)' 



^ ^exp 

a* 



mujr 



c-(l + m)ln(c) (l + o(l)) 
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T^u f f T Dci / ( f i + rn\ V , exp(-6HALL,d • c) , , 

Therefore for large c, PFI ^ — exp - c (1 + o(l)) = ■ (1 + oil)), 

a* V V 2 n J J ^HALL.d 

where a* = min (a ■ ITV.D^+I-^^^-^^^-^ 6hall d = i'muj./2) - (1 + m)/n, and 

^HALL.d = a*/K. 

Appendix C 

SADD FOR THE Boolean and the Path loss Models 

Fix 1 ^ i 5^ N. For each change time T ^ 1, define Tt = (^{xl''\s G A/", 1 ^ ^ T), and for 
£e G Ai, = (y{x\^\s eJ\fi,li^k^ T). From [9| (Theorem 3, Eqn. (24)), 

ess sup eJ.^^^^^)) (^(^™ie,{M) _ ^(1 + o(l)), as c ^ oo, (28) 

Define Ts^^,..k,{Mi)^T} the cr-field generated by the event {r™'^'^''^*) ^ T}, and similarly define the 
(7-field J|^ruiej,yj. Evidently J^|^ruie,(,)^y} ^ -^(^'^^_^-) and -^{r'-i'j^r} ^ -^(T-i)- By iterated conditional 
expectation, 

ECd(4)) (^(^ruie,(M) _ r)+|J-^,..^y}) ^ esssupE?(^=)) (^(^-ie,(M) _ T)+|7-(T_i)) (29) 
We can further assert that 

Using this observation with Eqn. |29] and Eqn. |28l we can write, as c ^ oo, 

E(d(^=)) (^(^rule,{M) _ T)+| ^ ^(1 + o(l)) (30) 

Finally, E^^^^^^ _ T)+\r^u^. ^ t) 1^,,..^^} = E?^'^" ((r-i^.(^-) - r)+|^{...^^}) I^,..^^}. 

We conclude, from|30l that, as c ^ oo, E^^^^^'>'> ((r™'^'(-^') - T)+|r^"'^ ^ T) ^ f (1 + o(l)). 
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